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Abstract 

The phase diagram of the coupled sine circle map lattice shows spatio-temporal intermittency 
of two distinct types: spatio-temporal intermittency of the directed percolation (DP) class, and 
spatial intermittency which does not belong to this class. These two types of behaviour are seen to 
be special cases of the spreading and non-spreading regimes seen in the system. In the spreading 
regime, each site can infect its neighbours permitting an initial disturbance to spread, whereas 
in the non-spreading regime no infection is possible. The two regimes are separated by a line 
which we call the infection line. The coupled map lattice can be mapped on to an equivalent 
cellular automaton which shows a transition from a probabilistic cellular automaton (PCA) to a 
deterministic cellular automaton (DCA) at the infection line. Thus the existence of the DP and 
non-DP universality classes in the same system is signalled by the PCA to DCA transition. We 
also discuss the dynamic origin of this transition. 
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The identification of tlie universality class of spatiotemporal intermittency l|] in spatially 
extended systems has been a long standing problem in the literature. Early conjectures 
argued that the transition to spatio-temporal intermittency is a second order phase transi- 
tion, and the transition falls in the same universality class as directed percolation j^. This 
conjecture has become the central issue in a long-standing debate [3|, |j, l5|, which is still not 
completely resolved. 

Studies of the coupled sine circle map lattice have thrown up a number of intriguing 
observations of relevance to this problem jo]. This system has regimes of spatio-temporal 
intermittency (STI) with critical exponents which fall in the same universality class as di- 
rected percolation (DP), as well as regimes of spatial intermittency (SI) which do not belong 
to the DP class. Both these regimes lie on the bifurcation boundaries of the synchronised 
solutions of the map. The spatio-temporally intermittent regime seen here has an absorb- 
ing laminar state, i.e. a laminar site remains laminar unless infected by a neighbouring 
turbulent site. The burst states spread and can percolate through the entire lattice. The 
system shows a convincing set of directed percolation exponents in this regime Gj. In the 
spatially intermittent regime, the laminar sites are frozen in time and the burst sites show 
temporally periodic or quasi-periodic behaviour. The laminar sites do not get infected by 
neighbouring turbulent sites. Hence, the spatially intermittent state is non-spreading and 
does not show directed percolation exponents. Thus, both DP and non-DP behaviour can 
be seen for different parameter regimes of the same system. 

In the present paper, we show that the infective directed percolation behaviour of spatio- 
temporal intermittency and the non-infective behaviour of spatial intermittency are special 
cases of the more general spreading to non-spreading transition seen in this system. The 
spreading and non-spreading regimes are separated by a line which we call the infection line. 
Above the infection line, the burst states can infect neighbouring laminar states and spread 
through the lattice, whereas below this line the burst states cannot infect their neighbours 
and the non-spreading regime is seen. The infection line intersects the bifurcation boundary 
of the synchronised solutions. Intermittent solutions are seen along this boundary, with the 
DP type of STI being seen above the infection line, and the non-DP SI being seen below 
the infection line. Spreading and non-spreading solutions are also seen off the bifurcation 
boundary. However, the distribution of laminar lengths shows power-law scaling only for 
parameter values which are very close to the bifurcation boundary, and falls off exponentially 
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FIG. 1: shows the phase diagram of the CML. Spatiotemporal intermittency of the DP class 
is seen at the O-s. Spatial intermittency with quasi-periodic bursts is seen at the A-s, SI with 
periodic bursts at the x-s and*; STI with travelling wave laminar states and solitons at the boxes 
(□). The spreading and the non-spreading regimes separated by the infection line are shown. The 
region of cross-over from STI of the DP class (indicated by O-s), to SI with quasi-periodic bursts 
(indicated by A-s) is shown in the inset. 



as the parameter values get more distant from the bifurcation boundary. Other exponents 
associated with DP behaviour are also observed only along the bifurcation boundaries. Fur- 
ther insights into the spreading to non-spreading transition are obtained by mapping the 
coupled map lattice (CML) onto a cellular automaton. The spreading to non-spreading 
transition seen across the infection line maps on to a transition from a probabilistic cellu- 
lar automaton to a deterministic cellular automaton. Thus the change from spreading to 
non-spreading behaviour seen in the CML is reflected in this transition. We also provide a 
pointer to the dynamic origin of this transition. 

The coupled sine circle map lattice studied here is known to model the mode-locking 
behaviour [Tj] seen in coupled oscillators, Josephson Junction arrays etc. The model is 



defined by the evolution equations 



= (1 - e)/(x*) + -ifixU) + f{xl,)] (mod 1) 



(1) 



where the index z is a discrete site index which runs on a one dimensional lattice of 
sites, and t s a discrete time index. The parameter e is the strength of the couphng between 



FIG. 2: The space time plots of (a) STI with synchronized laminar state and turbulent bursts 
seen at Q = 0.06, e = 0.7928. (b) SI with synchronized laminar state and quasi-periodic bursts seen 
at = 0.031, e = 0.42. 

the site i and its two nearest neighbours. The local on-site map, f{x) is the sine circle map 
defined as f{x) = x + Vl — ^ sin(27rx), where, K is the strength of the nonlinearity and fl 
is the winding number of the single sine circle map in the absence of the nonlinearity. We 
study the system with periodic boundary conditions in the parameter regime < Q < ^ 
(where the single circle map has temporal period 1 solutions), < e < 1 and K = 1.0. The 
phase diagram of this model evolved with random initial conditions is shown in Fig. [H 

The synchronised fixed point solutions = sin~^(^^) are seen in the regions indicated 
by dots in the phase diagram. The infection line seen in the figure separates the remaining 
part into the spreading and non-spreading regions. The space-time plots of the solutions 
seen in both these regions show co-existing laminar states and turbulent states (See Fig. 

The power law behaviour of the distribution of laminar lengths seen at the DP points 




FIG. 3: The space time plots seen (a) above the infection line at O = 0.062, e = 0.272 and (b) 
below the infection line at 17 = 0.062, e = 0.258. 
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FIG. 4: shows the laminar length distribution of (a) STI of the DP class seen at Q = 0.06, e = 
0.7928 and (b) SI with quasi-periodic bursts at 0, = 0.04, e = 0.402. The data have been obtained 
for a 10^ site lattice and are averaged over 50 initial conditions. 

(indicated by diamonds on the bifurcation boundary above the infection hne) is shown in 
Fig. |4](a). This scales with an exponent C = 1.67, characteristic of directed percolation. 
A full set of directed percolation exponents is seen at these points j^. The exponential 
fall off of the laminar length distributions for the spreading solutions off the bifurcation 
boundary, is shown in Fig. Hl^b). Below the infection line, in the non-spreading regime, 
the laminar sites are the synchronised fixed point x* and the burst sites can be temporally 
frozen, periodic or aperiodic. The power-law scaling of the laminar length distribution 
for spatial intermittency is also shown in Fig. H] (a). Here the distribution scales with 
exponent C = 1.1, distinct from the DP value. The exponential decrease of this distribution 
seen for other non-spreading solutions at points off the bifurcation boundary is shown in 
Fig. 111(b). Thus the sine circle map CML shows a transition from a spreading regime to 
a non-spreading regime. In order to gain further insights into this transition, we map the 
CML to a stochastic model, a probabilistic cellular automaton of the Domany-Kinzel type 







The equivalent cellular automaton, defined on a one dimensional lattice of size A^, is set 
up to mimic the dynamics of the laminar and burst states. The state variable f * at site i 
and at time t takes values vj = if the site is in the laminar state, and vj = 1 if the site is 
in the burst state. By the CML evolution equation ([1]), the state of the variable at site i at 
time t + 1 depends on the state of the variables at sites i, i — 1 and z + 1 at time t. Hence, 
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the probability of the site i at time t + 1 being in the burst state depends on the state of the 
sites i — and i + 1 at time t. We therefore define the CA dynamics in this system by the 
conditional probability F(v*'*"^|v*_i, , ^i+i). There are 2^ possible configurations and the 
symmetry between the sites i — 1 and z + 1 in the CML equation is used to obtain the effective 
probabilities p^'s, and define the CA rules as po = F(1|000) , = P(l|001) = F(l|100), 
P2 = P(l|010), P3 = P(l|011) = P(l|110), p4 = F(l|101) and p5 = 

The direct connection between equivalent CA and the CML of eq. ([1]) is set up by 
estimating the probabilities pk from the numerical evolution of the CML from random initial 
conditions from a uniform distribution over the (0, 1) interval for a given set of parameter 
values. The probabilities are estimated by finding the fraction of sites i which are in the 
burst state vj'^'^ = 1 at at time given that the site i and its nearest neighbours i — 1 and 
i + 1 existed in state k at time t. That is, the probability pk is estimated using pk = -tftttw^-i 
where and are the number of sites which at time t + 1 exist in the laminar states 
(f •"'""^ = 0) and the burst states (f*"*"^ = 1) respectively and at time t were the central sites 
of the configuration k. These probabilities, were extracted from a CML of size N = 2000 
averaged over 20000 timesteps discarding a transient of 30000 timesteps. The probabilities 
Pk in the spreading and non-spreading regimes in the phase diagram are listed in the Table 
[H It is clear that that the probability pq = P(1|000) is equal to zero in both regimes. This 
is the condition for an absorbing state where a laminar site with two laminar neighbours 





n 


e 


Po 


Pi 


P2 


P3 


P4 


Pb 


S(DP) 


0.060 
0.073 


0.7928 
0.4664 


0.0 
0.0 


0.220 
0.150 


0.0 
0.0 


0.933 
0.938 


0.627 
0.439 


0.984 
0.993 


S 


0.070 


0.264 
0.248 


0.0 
0.0 


0.140 
0.050 


0.0 
0.0 


0.982 
0.989 


0.391 
0.160 


0.999 
0.999 


NS 


0.070 


0.232 
0.228 


0.0 
0.0 


0.000 
0.000 


0.0 
0.0 


1.000 
1.000 


0.000 
0.000 


1.000 
1.000 


NS(SI) 


0.031 
0.044 


0.420 
0.373 


0.0 
0.0 


0.000 
0.000 


0.0 
0.0 


1.000 
1.000 


0.000 
0.000 


1.000 
1.000 



TABLE I: shows the probabilities pkS obtained in the spreading (S), non-spreading (NS) regimes 
and at DP and SI points averaged over 50 initial conditions. 
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cannot spontaneously evolve into a burst state. We also see that p2 = P(l|010) = so that 
a burst site with two laminar neighbours always goes into a laminar state, i.e. the laminar 
neighbours suppress the central burst site. The probabilities pi = P(l|001) = P(l|100) and 
P4 = P(l|101) are essentially infection probabilites by which a laminar site is infected by its 
burst neighbour or neighbours to change to a burst site. 

It is clear from Tabled that these probabilities show drastically different behaviour in the 
spreading and non-spreading regimes. In the case of the spatio-temporal intermittency with 
directed percolation exponent, which lies in the spreading regime of the phase diagram, these 
infection probabilities pi and p^ lie in the open interval (0, 1). Therefore, in the spreading 
regime, the dynamics is described by a probabilistic cellular automaton wherein the CA 
rules are probabilistic in nature. 

In the case of spatial intermittency which lies in the non-spreading regime, the probabili- 
ties obtained take the values or 1. In addition to po and p2 which are zero for the STI of the 
DP type, the infection probabilities pi and p^ go to zero, and p^ = P(l|011) = P(l|110) and 
P5 = P(l|lll) take the value 1. We also note that p2 is zero in the non-spreading regime, as a 
single burst site with two laminar neighbours is never observed in the non-spreading regime. 
Thus the probability of infection of a laminar state by its neighbouring burst state is zero in 
this regime, and no spreading of bursts can occur here. Cellular automata with probabilities 
which take values and 1 alone are called deterministic cellular automata (DCA), as given 
the state of the system at a time t, its state at the time t -|- 1 is deterministically known. 
Thus the spatial intermittency in the non-spreading regime can be represented by a DCA, 
up to the coarse graining defined earlier. Similar behaviour, either PCA (in the spreading 
regime) or DCA (in the non-spreading regime), is observed at other parameter points. 
A simple mean field approximation can be set up for the PCA [9]. Let rrit and rrit+i be 
the density of burst states in the lattice at the t^^ and t + 1*^ timestep. Using the CA 
rules defined above, the mean field evolution equation for the density of bursts is given by 
mt+i = {2pi +p2)mt{l - mtY + (2p3 + Pa)^^^^ - rrit) +P5m^. 

Approximating the values of Table [T] by p^ = 1, and using p2 = 0, the evolution equation 
reduces to m^+i = 2pi mf(l — rrit)^ + (2p3 -|- p^) mf{l — rrit) + nT't- This equation has 
three fixed points m = 0, m = ^^2pf-2p3-P4 ~ ~ The stability regions of these 

fixed points, as well as the co-existence region, where both the fixed points m = and 
m = 1 are stable, are shown in Figure O The DP probabilities seen in the spreading 
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FIG. 5: (Color online) shows the {pi,P3,P4} probabilities of the DP points (marked with diamonds) 
and for points in the spreading regime at = 0.065 (dotted line). The triangle corresponds to the 
probabilities in the non-spreading regime. 

regime, as well as the PCA probabilities seen at other points in the spreading regime, lie 
in this co-existence region with the 2ps + p4 = 2 plane as a lower bound. However, all the 
probabilities associated with the deterministic cellular automaton seen in the non-spreading 
regime (i.e. pi = 0, = 1, P4 = 0) lie at the vertex of the co-existence region in this cube 
(marked with a triangle in Fig. [5]). As the parameters Q and e vary along curves which 
cross the infection line, the PCA probabilities seen in the spreading regime collapse to the 
DCA probabilities at the vertex of the cube. Thus a PCA to DCA transition occurs at 
the infection line. Since the DP behaviour and the spatially intermittent behaviour along 
the bifurcation boundaries are special cases of spreading and non-spreading behaviour, the 
transition from the DP universality class of spatio-temporal intermittency to the non-DP 
universality class of spatial intermittency is reflected in the transition from the PCA to the 
DCA Q. 

The dynamical reason for this transition can be found by investigating the bifurcation 
diagram of the system (Fig. [6]). The range of e values on the vertical axis of Fig. [1] 
cut across the infection line at e = 0.254. The bifurcation diagram clearly shows that an 



attractor widening crisis 



ll| appears at this point. Similar behavior is seen for other sites. 



The spreading regime seen in the phase diagram emerges exactly at the point at which the 
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FIG. 6: shows the bifurcation diagram at 17 = 0.065 where x values for a typical site i have been 
plotted as a function of e, an initial transient of 14, 000 iterates has been discarded. 

attractor widens, with the non-spreading regime corresponding to the pre-widening regime. 
This widening also identifies the point at which the equivalent cellular automaton undergoes 
a PCA to DCA transition. In the pre-crisis region, each site follows either a periodic or 
quasiperiodic trajectory and is not infected by the behaviour of its neighbours. Thus, its 
CA analogue is deterministic as listed in Table [B In the postcrisis regime, each site is able 
to access the full x range, as well as infect its neigbours,and the bursting and spreading 
behaviour characteristic of the spreading regime is seen. This is reflected in the equivalent 
cellular automaton by a transition to PCA behavior (Table [I]). It is to be noted that the 
volume of the attractor in phase space will be much larger post-crisis, as compared to the 
pre-crisis volume. Further characterisation of the crisis is in progress. 

To conclude, the spatio-temporal intermittency of the directed percolation class and 
the spatial intermittency of the non-directed percolation class, seen along the bifurcation 
boundaries of synchronised solutions, are special cases of the spreading and non-spreading 
regimes seen off the bifurcation boundaries. The two regimes are separated by the infection 
line which intersects with the bifurcation boundary of the synchronised solutions at the 
point where the cross-over between the directed percolation and non-directed percolation 
behaviour takes place. Thus the behaviour seen in coupled sine circle map lattice is organised 
around the locations of the bifurcation boundaries and the synchronised solutions, and the 
infection line. The existence of two distinct universality classes, in the phase diagram of the 
sine circle map lattice is a reflection of the transition of the equivalent cellular automaton 
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from the probabilistic phase to the deterministic phase and the concomitant suppression of 
the spreading or infectious modes. The dynamic origins of this transition he in an attractor- 
widening crisis. We beheve this is the first time that such a direct connection has been found 
between a dynamical phenomenon viz. a crisis in an extended system and the statistical 
properties of the extended system viz. the exponents and universality classes. Similar 
directed percolation to non-directed percolation transitions have been seen in other coupled 
map lattices, as well as in pair contact processes, solid on solid models and models of non- 
equilibrium wetting 12|. Our results may have useful pointers for the analysis of other 
systems, and thus contribute to the on-going debate on the identification of the universality 
classes of spatiotemporal systems. 
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